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Abstract 

Annihilation of dark matter particles accumulated in the Sun would produce a flux of high-energy 
neutrinos whose prospects of detection in neutrino telescopes and detectors have been extensively 
discussed in the literature. However, for many channels into Standard Model particles, there would 
also be a flux of neutrinos in the MeV range from the decays at rest of muons and positively charged 
pious. These low-energy neutrinos have never been considered before and they open the possibility 
to also constrain dark matter annihilation in the Sun into e~^e^ , fJ'^fJ^^ or light quarks. Here we 
perform a detailed analysis using the recent Super-Kamiokande data in the few tens of MeV range 
to set limits on the WIMP-nucleon scattering cross section for different annihilation channels. 
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Introduction. — 

There is overwlielming evidence of tlie existence of a massive non-baryonic dark compo- 
nent which contributes to about 80% of the energy budget of the Universe being a weakly 
interacting massive particle (WIMP), with mass lying from the GeV to the TeV scale, one 
of the most popular candidates. 

One of the different proposed strategies to detect WIMPs is to search for the flux of 
high-energy neutrinos from the annihilations of WIMPs accumulated in the center of the 
Sun |2|. Many different studies have^valuated the prospects of detection of these neutrinos 
with neutrino telescopes/detectors [3]. However, previous works have focused on WIMPs 
annihilations into hadronic or r+r~ channels. On the other hand, annihilations into 
or light quarks have always been neglected as muons and pions lose energy very effectively in 
the dense regions where they would be produced and then would decay at rest, giving rise to 
neutrinos in the MeV range. Likewise, annihilations into e~^e~ have never been considered, 
as e^/e~ would not produce directly neutrinos. Nevertheless, in their propagation they 
could interact with nuclei and produce pions, which would be stopped. The tt^ would then 
get captured and subsequently absorbed by the nuclei of the medium without decaying, 
but the 7r+ would decay at rest, producing as well a flux of MeV neutrinos. On the other 
hand, hadronic and t^t" channels, along with heavy mesons (the source of the high-energy 
neutrinos considered so far), would also produce pions. The energies of these neutrinos lie 
at the energy range where the diffuse supernova neutrino background (DSNB) is searched 
for by detectors as Super-Kamiokande (SK). 

Here we consider, for the first time, the potential signal of these low-energy neutrinos 
from WIMPs annihilations in the Sun and use the most recent SK data [4] to set bounds 
on the scattering cross section of WIMPs off nucleons for different annihilation channels. 



MeV neutrinos from WIMPs annihilations in the Sun. — 

Galactic WIMPs could get gravitationally bound within the Sun if after many elastic 
scatterings off solar nuclei they lose energy and their velocity gets lower than the Sun's 
escape velocity. Typically, WIMPs would make many orbits between interactions, so their 
density could be approximated as an isothermal sphere, n^{r) oc e~'^^'^^'^'>^'^^ Pj^I, with 0(r) 
the gravitational potential and the WIMPs temperature, calculated following Ref. j5| and 
using the Standard Solar Model (SSM) ^. 
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The evolution of the number of WIMPs in the Sun depends on the capture rate, C©, on 
the annihilation cross section times the relative velocity per unit volume, Aq, and on the 
evaporation rate, Eq, which is only relevant for low- mass WIMPs. We use the capture rates 
calculated in DarkSUSY Q| and we calculate Aq and Eq following Ref. [9i] and using the 
SSM [7]. The WIMPs annihilation rate now {t = Iq = 4.57 Gyr) is given by jol 

1 f tanh(/tt0/rg) ^ 
r(m„a,)--Ce + i tanh(«: t/r^) J ' 

where te = {AqCq)~^^'^ is the equilibration time scale in the absence of evaporation and 

we assume an annihilation cross section typical 
of thermal WIMPs, {a^v) = 3 ■ lO'^^cm^/s and, for the scattering cross sections under 
consideration, equilibrium is always reached (tQ » te, tanh(fi; tQ/r^;) ^ 1). 

Being produced in a very dense medium, among all the final products of WIMPs annihila- 
tions, only neutrinos can escape. So far, all previous works have focused on the high-energy 
neutrino flux resulting from the subsequent hadronization, fragmentation and decay of the 
final states in heavy quarks, gauge bosons or r+r^ channels, and have disregarded annihi- 
lations into e~^e~, ^~ or light quarks because they would only produce (if any) a flux of 
low-energy neutrinos from pion and muon decay at rest. These MeV neutrinos are the focus 
of this work. 

The propagation of the annihilation products would produce pions that would be stopped 
and could subsequently decay at rest, giving rise to a monochromatic neutrino spectrum at 
29.8 MeV. Practically all vr", after stopping, are captured in an atomic orbit and get absorbed 
by the nucleus without decaying. Hence, only neutrinos from decays would contribute. 
Let us note that the propagation of high-energy e~ /e^ would also produce small amounts 
of pions. In addition, all muons produced in pion decays, in the leptonic decay modes of 
quarks or taus and in the case of direct annihilations into fi~^fJ^~, are stopped and decay at 
rest. Thus, in addition to the monochromatic spectrum from 7r+ decays at rest, neutrinos 
from /i+ and fi~ decaying at rest would also contribute (with a well known spectrum below 
52.8 MeV) to the final low-energy neutrino flux. 

In order to calculate the relative contributions of each type of neutrino spectrum, we 
simulate all the particles propagation with GEANT4 jlO|. To determine for each WIMP 
mass the average density and composition of the medium where the products of WIMPs 
annihilation propagate, we take the SSM Q and use the WIMP distribution in the Sun. 



For the simulations, we proceed as follows. For the case of WIMPs annihilating into a 
pair of leptons we inject the two leptons with energies equal to the WIMP mass directly 
into GEANT4 and let them propagate and decay. In the case of WIMPs annihilations into 



quarks, we use PYTHIA 6.4 ll| to hadronize and fragment the initial quarks and do not let 
decay any of the final particles that are produced. Then, we inject into GEANT4 the full 
spectrum of all the produced particles and simulate their propagation in the Sun. Finally, 
we count the number of 7r+, tt", /i+ and fi~ that decay at rest, which allow us to compute 
the initial electron and muon neutrino and antineutrino fluxes at the source. 

These neutrinos would then propagate from the Sun to the Earth and would be detected 
via the charged- current interactions of those arriving at the detector in the electron flavor. 
In order to calculate the electron neutrino and antineutrino fluxes at the Earth, one should 
take into account that neutrinos mix. For the energies of interest (above 10 MeV) neutrinos 
propagate adiabatically in the Sun and at the region where they are produced matter effects 
are dominant, so (z/g) exit the Sun as almost purely z/2 (z/i) and z/^ (z/^) almost as an equal 
mixture of ui {U2) and u^. Hence, the probabilities for neutrinos to arrive at the Earth in 
the electron flavor are 

P{^^J. ^e) = cos^ 6*23 cos^ 6^12 cos^ 6113 + sin^ cos^ sin^ ^12 (1 + sin^ 612) 

P{^fi ~^ ^e) = cos^ 6*23 sin^ 9i2 cos^ 9i3 + sin^ 613 cos^ 9i3 sin^ 612 (1 + cos^ 612) 

P{h'e z^e) = sin^ 9i2 cos^ 6'i3 + sin^ ^13 

P{Ve J^e) = COS^ 6*12 COS^ 6*13 + slu"^ 6*13 . (2) 

The final electron neutrino and antineutrino spectra at the detector are obtained by 
combining the initial fluxes with the above probabilities (and we have neglected the small 
correction due to the Earth- matter effect). We use the values of the mixing angles from 
Ref. y. 

Detection of MeV neutrinos with SK. — 

SK is a water Cerenkov detector with a fiducial volume of 22.5 ktons. For energies below 
52.83 MeV, the best present data come from the search for the DSNB which is split into 
three phases: SK-I (tj = 1497 days), SK-II {tn = 794 days) and SK-III {t^J^ = 562 days). 
The recent analysis 1^ has substantially improved over the previous one 
efficiencies, lower energy thresholds and almost twice as much statistics. 
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The signal at SK is the detection of the positrons (electrons) produced in z/g (i^e) charged- 
current interactions in the detector below ~ 100 MeV. At these energies, the inverse beta- 
decay cross section {Uf,p — ne~^) is by two orders of magnitude larger that the u — e elastic 
scattering cross section. Below ~80 MeV, this is the dominant interaction of z/e- However, 
the interactions z/g and z/g off oxygen nuclei give non-negligible contributions and need also 
to be taken into account. 

The low energy threshold in the analysis is determined by the ability to remove the 
radioactive spallation caused by cosmic-ray muons hitting an oxygen nucleus. An improved 
spallation cut has allowed to reduce the energy threshold used in the previous analysis down 
to 16 MeV (17.5 MeV) for SK-I/III (SK-II). In addition to this, a number of other cuts were 
performed, as noise reduction, fiducial volume, solar angle, incoming event, decay electron, 
pion, Cerenkov angle and other cuts ^. The maximum energy in the recent SK analysis is 
88 MeV, which is also higher than in the previous one. 

In this energy range, the two dominant backgrounds which remain after the cuts are 
the atmospheric z/g and z/g background and, mainly, the Michel positrons (electrons) from 
the decays of low energy muons, produced by atmospheric z/^ (z/^) with typical energies 
of about ~200 MeV, which are below detection threshold, the so-called invisible muons. 
These muons are slowed down rapidly and subsequently decay, mimicking the signal from 
z/g or Uf., but with a spectrum whose shape is very well known. 

Analysis. — 

n 

In this work, we follow closely the analysis performed in Refs. although updating 
some aspects following the new SK analysis j^]. 

We have considered both the interactions of z/g off free protons and the interactions of 
z/g and ^off bound nucleons. For the inverse beta-decay reaction we use the full cross 
section |l5| and for the interactions off bound nucleons, we consider a relativistic Fermi gas 
model 16| with a Fermi surface momentum of 225 MeV and a binding energy of 27 MeV. For 
each of the three SK phases, we have used the corresponding energy-dependent efficiencies ^ 
and a Gaussian energy resolution function \v\. 

In order to obtain the upper limit on the WIMP-nucleon scattering cross section, we 

n 

use the recent data reported by the SK collaboration [4]. Unlike the SK analysis, we only 
consider the two backgrounds described above and only take the data in the Cerenkov angle 
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region defined as the signal region We expect this to affect very httle our results because 
the two extra sources of background considered in the SK analysis are subdominant at 
the percent level and their energy distribution is very flat, while the normalizations of the 
potential signal and of each background are left free in our analysis, which are the three 
parameters we fit (a, /3 and 7, which represent the total number of signal and background 
events). We consider 18 4-MeV bins (in the interval 16-88 MeV) and perform an unbinned 
maximum likelihood fit. As in the SK analysis, we obtain the best fit as the combination of 
parameters that maximizes the likelihood, defined as 

18 

£ = -Q [(« . Ai) + (/3 ■ Bi) + (7 ■ Ci)f^ , (3) 

1=1 

where the sum I is over all energy bins, Ni is the number of detected events in the /th 
bin, and A/, Bi and Ci are the fractions (so that for each case the total is normalized to 
1) of the WIMPs annihilation signal, Michel positrons and electrons from muon decay and 
atmospheric and Vf. spectra that are in the /th bin, respectively. The fractions Ai are 
calculated using the z/g and z/g low-energy fluxes from WIMP annihilations in the Sun as 
described above. When computing Bi we have taken into account that 18.4% of the /i~ are 
stopped in water and the spectrum of their decay electron gets distorted. In the calculation 
of the fractions Bi and Ci we have used the atmospheric neutrino flux calculation with 
FLUKA [18|. 

We have also included the energy-independent efficiency systematic error by modifying 
the likelihood in the way described in Ref. Q|, with a total error different for each of the 
data-taking phases. The final likelihood is maximized for each SK phase separately, so that 
it remains as a function of just a, the number of signal events. Finally, the three likelihoods 
are calculated as a function of the number of signal events/year and multiplied. The 90% 
confidence level (CL) limit on the number of signal events/year is determined by 

r°^° £tot(tt) da , , 

:iO_ ' = 0.9 . (4) 

/o°° Aot(a) 

The limit on the scattering cross section (cr^°) is then obtained by solving 



r(m„af)A- = a9o , ^ , (5) 



where A^^ is the number of events per WIMPs annihilation for each SK phase at the 
detector. 
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Figure 1: Limits on the spin- dependent scattering cross section of WIMPs off protons 
at 90% CL for different annihilation channels. The limits from SIMPLE, PICASSO and 



CDMS 



19j are shown with black lines. The limits from SK searches of GeV neutrinos are depicted 



for two annihilation channels {20/. The DAMA/LIBRA 123] regions (at 90% CL and 3a CL) are 



0/ 



also shown as interpreted in Refs. 



The results for spin-dependent (off protons) and spin-independent cross sections are 
shown in Figs. [1] and 121 respectively, along with the bounds from direct detection searches 
and the analysis of GeV neutrinos at SK. The most stringent bounds we obtain are for 
WIMPs annihilations into light quarks, which is the channel that generates the largest 
amount of pions. Whereas for spin- independent, the limits are one or two orders of 
magnitude weaker than the ones obtained with direct searches, for spin-dependent and low 
WIMP mass, they are comparable to (or even more stringent than) them. We note that, in 
this case, the limits for WIMPs annihilations into quarks exclude the DAMA/LIBRA |2l| 
region at 90% CL. We also obtain bounds for the case of WIMPs annihilations into e^e~ , 
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Figure 2: Limits on the spin-independent scattering cross section of WIMPs at 90% CL 
for different annihilation channels. The limits from CDMS, COGENT, XENON 10 and 
XENON 100 /23| / are shown with black lines. The limits from SK searches of GeV neutrinos are 



depicted for two annihilation channels 



'2Q]. The DAMA/LIBRA IM] regions (at 90% CL and 



3(7 CL) are also shown as interpreted in Refs. 



22] 



although very weak. The hmit for the jj,^ ji^ channel follows the dependence of the capture 
rate with the WIMP mass, as the initial muons are always stopped. On the other hand, for 
the r+r~ channel, at low masses, the number of muons decaying at rest is smaller because 
a fraction of taus goes to electrons and all the produced tt^ get absorbed without decaying. 
However, the passage through the solar medium of the products of tau decay gives rise to 
a number of pions that increases with energy and hence, with the WIMP mass (the same 
trend is seen for the case of WIMPs annihilations into quarks). This explains why these 
limits cross each other. 
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Conclusions. — 

The potential signal of GeV neutrinos produced after WIMPs annihilations in the Sun 
has been extensively studied so far [3] . However, these searches do not consider annihilations 
into light quarks, ^'^ ^~ or e^e~ , as these channels do not produce GeV neutrinos. However, 
the propagation in the Sun of the products of WIMPs annihilations into quarks or leptons 
would always produce pions and muons that after being stopped would decay at rest. In 
this work, we have considered these MeV neutrinos from pion and muon decay at rest, 
which represents a novel way of constraining WIMPs annihilations in the Sun. In order 
to do so, we have used the SK data on the DSNB ^ and have performed and analogous 
analysis. Our results. Figs. [UandEl show that these new limits are complementary to those 
from direct searches, mainly for low-mass WIMPs. 



Note added: As this work was being completed, we learned of an i ndep endent study 
with similar results by C. Rott, J. Siegal-Gaskins and J. F. Beacom |24| that will be 
submitted to arXiv simultaneously. 

Acknowledgment s . — 

We thank K. Bays and C. Pena-Garay for providing us with useful information, 
M. Vicente- Vacas for discussions and M. Sorel for encouragement. SPR thanks the Galileo 
Galilei Institute for Theoretical Physics, where parts of this work were done, for its hos- 
pitality. NB is supported by the DFG TRR33 'The Dark Universe'. JMA is supported 
by the Spanish Grant CONSOLIDER-Ingenio 2010 CSD2008-0037 (CUP) of the MINECO. 
SPR is partially supported by the Portuguese FCT through CERN/FP/123580/2011 and 
CFTP-FCT UNIT 777, which are partially funded through POCTI (FEDER) and by the 
Spanish Grant FPA2011-23596 of the MINECO. 



[1] G. Bertone, D. Hooper and J. Silk, Phys. Kept. 405, 279 (2005) [arXiv:hep-ph/0404175]. 
[2] J. Silk, K. A. Olive and M. Srednicki, Phys. Rev. Lett. 55, 257 (1985); M. Srednicki, K. A. Olive 
and J. Silk, Nucl. Phys. B 279, 804 (1987); T. K. Gaisser, G. Steigman and S. Tilav, Phys. 



9 



Rev. D 34, 2206 (1986). 

[3] M. Kamionkowski, Phys. Rev. D 44, 3021 (1991); A. Bottino et al, Phys. Lett. B 265, 

57 (1991); F. Halzen, T. Stelzer and M. Kamionkowski, Phys. Rev. D 45, 4439 (1992); 

L. Bergstrom, J. Edsjo and P. Gondolo, Phys. Rev. D 55, 1765 (1997) [hep-ph/9607237]; 

V. D. Barger et al, Phys. Rev. D 65, 075022 (2002) [hep-ph/0105182]; M. CirelH et al, and 

F. Vissani, Nucl. Phys. B 727, 99 (2005) [Erratum-ibid. B 790, 338 (2008)] [hep-ph/0506298]; 

O. Mena, S. Palomares-Ruiz and S. PascoH, Phys. Lett. B 664, 92 (2008) [arXiv:0706.3909 

[hep-ph]]; C. R. Das et al, arXiv:1110.5095 [hep-ph]. 
[4] K. Bays et al. [Supcr-Kamiokande Collaboration], Phys. Rev. D 85, 052007 (2012) 

[arXiv:1111.5031 [hep-ex]]. 
[5] D. N. Spergel and W. H. Press, Astrophys. J. 294, 663 (1985). 
[6] J. Faulkner and R. L. Gilliland, Astrophys. J. 299, 994 (1985). 

[7] A. M. Serenelh, W. C. Haxton and C. Pena-Garay, Astrophys. J. 743, 24 (2011) 

[arXiv: 1 104. 1639 [astro-ph.SR]] . 
[8] P. Gondolo et al, JCAP 0407, 008 (2004) [astro-ph/0406204] . 

[9] K. Griest and D. Seckel, Nucl. Phys. B 283, 681 (1987) [Erratum-ibid. B 296, 1034 (1988)]. 
[10] S. Agostinelh et al. [GEANT4 Collaboration], Nucl. lustrum. Meth. A 506, 250 (2003); J. Al- 
lison et al, IEEE Trans. Nucl. Sci. 53, 270 (2006). 
[11] T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605, 026 (2006) [hep-ph/0603175]. 
[12] D. V. Forero, M. Tortola and J. W. F. Valle, arXiv: 1205.4018 [hep-ph]. 

[13] M. Malek et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 90, 061101 (2003) [hep- 
ex/0209028]. 

[14] S. Palomares-Ruiz and S. Pascoli, Phys. Rev. D 77, 025025 (2008) [arXiv:0710.5420 [astro-ph]]; 

S. Palomares-Ruiz, Phys. Lett. B 665, 50 (2008) [arXiv:0712.1937 [astro-ph]]. 
[15] A. Strumia and F. Vissani, Phys. Lett. B 564, 42 (2003) [astro-ph/0302055]; P. Vogel and 

J. F. Beacom, Phys. Rev. D 60, 053003 (1999) [hep-ph/9903554]. 
[16] R. A. Smith and E. J. Moniz, Nucl. Phys. B 43, 605 (1972) [Erratum-ibid. B 101, 547 (1975)]. 
[17] J. Hosaka et al. [Super-Kamiokande Collaboration], Phys. Rev. D 73, 112001 (2006) [hep- 

ex/0508053]; J. P. Cravens et al. [Super-Kamiokande Collaboration], Phys. Rev. D 78, 032002 

(2008) [arXiv:0803.4312 [hep-ex]]; K. Abe et al. [Super-Kamiokande Collaboration], Phys. Rev. 

D 83, 052010 (2011) [arXiv:1010.0118 [hep-ex]]. 



10 



[18] G. Battistoni et al, Astropart. Phys. 19, 269 (2003) [Erratum-ibid. 19, 291 (2003)] [hep- 

ph/0207035]; and Astropart. Phys. 23, 526 (2005). 
[19] M. Felizardo et al. [SIMPLE Collaboration], Phys. Rev. Lett. 108, 201302 (2012) 

[arXiv: 1106.3014 [astro-ph.CO]]; S. Archambault et al. [PICASSO Collaboration], Phys. Lett. 

B 711, 153 (2012) [arXiv: 1202. 1240 [hep-ex]]; E. Behnke et al. [COUPP Collaboration], 

arXiv: 1204.3094 [astro-ph.CO]. 
[20] T. Tanaka et al. [Super-Kamiokande Collaboration], Astrophys. J. 742, 78 (2011) 

[arXiv: 1108.3384 [astro-ph.HE]]; R. Kappl and M. W. Winkler, Nucl. Phys. B 850, 505 (2011) 

[arXiv: 1104.0679 [hep-ph]]. 
[21] R. Bernabei et al. [DAMA/LIBRA Collaboration], Eur. Phys. J. C 56, 333 (2008) 

[arXiv:0804.2741 [astro-ph]]; and Eur. Phys. J. C 67, 39 (2010) [arXiv: 1002. 1028 [astro- 

ph.GA]]. 

[22] C. Savage et al, JCAP 0904, 010 (2009) [arXiv:0808.3607 [astro-ph]]; and Phys. Rev. D 83, 
055002 (2011) [arXiv: 1006.0972 [astro-ph.CO]]. 

[23] C. E. Aalseth et al. [CoGeNT Collaboration], Phys. Rev. Lett. 106, 131301 (2011) 
[arXiv:1002.4703 [astro-ph.CO]]; D. S. Akerib et al. [CDMS Collaboration], Phys. Rev. D 
82, 122004 (2010) [arXiv: 1010.4290 [astro-ph.CO]]; J. Angle et al. [XENONIO Collabora- 
tion], Phys. Rev. Lett. 107, 051301 (2011) [arXiv:1104.3088 [astro-ph.CO]]; E. Aprile et al. 
[XENONIOO Collaboration], arXiv: 1207.5988 [astro-ph.CO]. 

[24] C. Rott, J. Siegal-Gaskins and J. F. Beacom, submitted to arXiv. 



11 



